DECAY OF THE VLASOV-POISSON-BOLTZMANN SYSTEM 
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Abstract. We establish the time decay rates of the solutions to the Cauchy problem for the 
two-species Vlasov-Poisson-Boltzmann system near Maxwellians via a refined pure energy method. 
The negative Sobolev norms are shown to be preserved along time evolution and enhance the 
«) ■ decay rates. The total density of two species of particles decays at the optimal algebraic rate as 

' the Boltzmann equation, but the disparity between two species and the electric field decay at an 

exponential rate. This phenomenon reveals the essential difference when compared to the one- 
species Vlasov-Poisson-Boltzmann system in which the electric field decays at the optimal algebraic 
rate or the Vlasov-Boltzmann system in which the disparity between two species decays at the 
optimal algebraic rate. Our proof is based on a family of scaled energy estimates with minimum 
QQ ' derivative counts and interpolations among them without linear decay analysis, and a reformulation 

^SJ , of the problem which well displays the cancelation property of the two-species system. 



> 
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1. Introduction 



< 

"j^ I The dynamics of charged dilute particles (e.g., electrons and ions) in the absence of magnetic 

effects can be described by the Vlasov-Poisson-Boltzmann system: 

dtF+ + V ■ V,F+ + • V,F+ = Q(F+, F+) + Q(F+, F_), 

^, 5tF_+^;-V,F_ -V,$-V,F_ = Q(F_,F+) + Q(F_,F_), 

in 

m 
cn 

I with initial data F-i-{0,x,v) = Fq^±{x,v). Here F±(t,x,v) > are the number density functions 

for the ions (-|-) and electrons (— ) respectively, at time t > 0, position x = (xi,X2,X3) G M'^ and 
velocity v = {vi,V2,V3) £ M^. The self-consistent electric potential <I> = ^{t,x) is coupled with 
F±{t,x,v) through the Poisson equation. The collision between particles is given by the standard 
Boltzmann collision operator Q{hi,h2) with hard-sphere interaction: 



A,.$= / F<-F_dv, 



Q{hiM){v)= / \{u - v) ■ u\{hi{v')h2{u') - hi{v)h2{u)} du. (1.2) 
Here w G and 

v' = V — [{v — u) ■ wjo;, u' = u+[{v — u) ■ wjo;, (1.3) 

which denote velocities after a collision of particles having velocities v and u before the colli- 
sion, and vice versa. Since the presence of all the physical constants does not create essential 
mathematical difficulties, for notational simplicity, we have normalized all constants in the Vlasov- 
Poisson-Boltzmann system to be one. Accordingly, we normalize the global Maxwellian as (with 
= 0) 

^,{v) = ^+{v) = f,^{v) = ^-L_e-l^l'/2. (1.4) 
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In this paper, it is more convenient to consider the sum and difference of and F_. This is 
motivated by some previous works [H \TT\ [22] on the other systems of binary fluids. Defining 

F = F+ + F_ and G = F+-F^, (1.5) 

that is, F is the total density for the two species of particles and G represents the disparity between 
two species. Then the system (jl.ip can be written as the equivalent form: 

dtF + v V,F + V,^ ■ V,G = Q(F, F), 
dtG + vV,G + V,^-V,F = Q{G,F), 

A,$= / Gdv, 



with initial data F{0,x,v) = Fq{x,v) and G{0,x,v) = Gq{x,v). We define the standard perturba- 
tion [f{t,x,v),g{t,x,v)] around the corresponding equilibrium state [^,0] as 

F = fi + and G = ^/Jig, (1.7) 

then the Vlasov-Poisson-Boltzmann system for the perturbation [/, g] takes the form 

dtf + v V,./ + £i/ = 9Ti := r(/, /) + iv,.cl> • vg - V,<I> • V,g, 

dtg + v- V^g - V^^> • + £.29 = ^2 := r{g, f) + ^V^^ • vf - • V„/, (1.8) 



A:,.^ = / g^dv, 



with initial data f{0,x,v) = fo{x,v) and g{0,x,v) = gQ{x,v). Here D^Ii and ^2 represent the 
nonlinear terms. The well-known linearized collision operator Ci and another linearized operator 
£2 are defined as 

£ih = --^{Q{n,^h) + Q{^h,fi)}, C2h = --^Q{^h,fi), (1.9) 

and the nonlinear collision operator (non-symmetric) is given by 

T{hi,h2) = -^Qi^hi, ^h2). (1.10) 

Note that the linear homogeneous system of (jl.Sp is decoupled into two independent subsystems. 
One is the Boltzmann equation for /, and the other one is a system almost like the Vlasov-Poisson- 
Boltzmann for g and but with a different linearized collision operator €2- 

Notice that [Ci , C2] is equivalent to the linearized collision operator L defined in [6] . It is well- 
known that the operators £1 and £2 are non-negative. For any fixed {t,x), the null spaces of £1 
and £2 are given by, respectively, 

7V(£i) = span{.y//i, VyZ/X; ^^"^ 7\A(£2) = spanj^Z/i}. (l-H) 

For any fixed {t,x), we define Pi as the orthogonal projection on the null space AA(£i). Thus 
for any function f(t,x,v) we can decompose 

/ = Pi/ + {I-Pi}/, (1.12) 

where Pi/ is called the hydrodynamic part of / and {I — Pi}/ is the microscopic part. We can 
further denote 

Pi/ = {ait, x) + b{t, x)-v + c{t, x) f ^ - ^/a^- (1-13) 
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Here the hydrodynamic field of /, [a{t,x),b{t,x),c{t,x)], represents the density, velocity and tem- 
perature fluctuations physically. Similarly, for any fixed (t, x), we define P2 as the orthogonal 
projection on the null space M{C2)- Thus for any function g{t,x,v) we can decompose 

g = F29 + {1-^2)9, (1.14) 

and we can further denote 

P2g = dit,x)^. (1.15) 
Here d{t, x) represents the concentration difference fluctuation. 

Notation. In this paper, with an integer ^ > stands for the usual any spatial derivatives 
of order £. When £ < or £ is not a positive integer, stands for defined by ()A.7p . We use 
H^(W^), s G M to denote the homogeneous Sobolev spaces on with norm ||-|| defined by (jA.Sp . 
and we use H^{W^) to denote the usual Sobolev spaces with norm \\-\\ffs and L^(M'^), 1 < p < 00 to 
denote the usual spaces with norm IHI^^p. 

We shall use (•, •) to denote the inner product in Rj^ with corresponding norm | • I2, while 
we use (•, •) to denote the inner product either in M.^ x or in with norm || • || without 
any ambiguity. We shall simply use L^, to denote L^(M^) and L^(M^) respectively, etc. We will 
use the notation L^H^ to denote the space L^(M^; H^) with norm 




(1.16) 



and similarly we use the notations of L^H^, L^Lx and LxL"^, etc. We define the space- velocity 
mixed derivatives by 

«7 _ P1I2 pil3 piPi 

'-'13 ~ ^Xi'-'x2'-'x3^Vi '~'V2 

where 7 = [71,72,73] is related to the space derivatives, while /3 = [/3i, /32, /^s] is related to the 
velocity derivatives. 

For the Boltzmann operator (II. 2p . we define the collision frequency as 

u^v) = \v — u\fi{u)du, (1-17) 

which behaves like 1 + We define the weighted norms 

I |2 I 1/2 |2 II ||2 II 1/2 ||2 /I io\ 

lalu = W ' 5I2. IbIL = w ' 9\\ ■ (1-18) 

We denote by the weighted space with norm || • 

Throughout this paper, C > will denote a generic constant that can depend on the parameters 
coming from the problem, and the indexes and s coming from the regularity on the data, but 
does not depend on the size of the data, etc. We refer to such constants as "universal." Such 
constants are allowed to change from line to line. We will employ the notation a < 6 to mean 
that a < Cb for a universal constant C > 0. We also use Co for a positive constant depending 
additionally on the initial data. To indicate some constants in some places so that they can be 
referred to later, we will denote them in particular by Ci,C2, etc. 

As it can be obtained by formally setting the magnetic field to be zero in the two-species Vlasov- 
Maxwell-Boltzmann, the global existence of classical solutions near Maxwellians to the two-species 
Vlasov-Poisson-Boltzmann system (II. ip can be found in [6l [Tl] for the spatially periodic domain 
and in [T9] for the whole space. For the one-species Vlasov-Poisson-Boltzmann system, the first 
global unique solution near Maxwellians was constructed in [5] for the periodic domain and in j25] 
for the whole space under the restrictions that either the mean free path is sufficiently small or the 
constant background charge density is sufficiently large. Those restrictions in [25] were removed 
in [26], and [4] proved the global existence of solutions for more general nonconstant background 
density. These proofs are based on the nonlinear energy method developed in jSj [U [71 E] and in 
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[mni]. It seems that those results on the existence of solutions can be directly generalized to the 
two-species Vlasov-Poisson-Boltzmann system without any additional essential difficulties. 

The time decay rate of the solutions has been an important problem. It is well known in 
[6] that for the periodic domain the solutions to the Vlasov-Poisson-Boltzmann system decay at 
an exponential rate similarly as the Boltzmann equation [2QI [T7]. For the whole space, |26] also 
obtained a time decay rate of 0(t~^/^) in the L^L^-norm of the solution for the one-species Vlasov- 
Poisson-Boltzmann system based on the pure energy method and a time differential inequality, and 
this result was generalized to the two-species system in |27j . It is noticed that the decay rates in 
[26] and [27] are not optimal. Concerning the optimal decay rates, under the additional assumption 
that L^L;^-norm of the initial perturbation /q is sufficiently small, by combining the linear decay 
results from the Fourier analysis and the nonlinear energy estimates, proved that the L^-norm 
of solutions to the one-species Vlasov-Poisson-Boltzmann system decay at the optimal rate of 
0(t~^/^) which is slower than the 0{t~^^^) of the Boltzmann equation. While under the additional 
assumption that L^L;^-norm of the initial perturbation /q = [/o,+ ,/o,-] and the L^-norm of Vx^o 
are sufficiently small, combining the construction of the compensating functions for two-species 
system by using Kawashima's method [12j and the nonlinear energy estimates, pl[j proved that the 
L^-norm of solutions to the two-species Vlasov-Poisson-Boltzmann system decay at the optimal 
rate of 0{t~^/^) which is same as that of Boltzmann equation \16\ 121] . One may also find from 
[3] that the L^L^ -norm with 2 < r < oo of solutions to the two-species Vlasov-Poisson-Boltzmann 
system decay at the optimal rate of 0(t~^/^*^^~^/'~^). We remark that if given the same condition 
on the initial data (the difference between [2j and |24^ [3] is whether one imposes the condition that 
||Va;$||/,i is small, which is equivalent to that || Va;A~^(/o, y^) is small), [2] and [211 [3] should 
get the same optimal decay rate of the solutions. It is observed that the electric field is coupled 
with / through the Poisson equation, so one may obtain the algebraic decay rate of V^^ from that 
of / in [2l|3l[23j. One may argue by contradiction that from the structure of the problem, the decay 
rate of electric field of the one-species system is optimal. However, we will show in this paper that 
for the two-species system the electric field decays at an exponential rate. This is totally due to 
the special cancelation property of the two-species system. 

It is difficult to show that the norm of the solution can be preserved along time evolution in 
the approach [H O [SI [Ml [H] . On the other hand, the existing pure energy method [26[ [27] 

of proving the decay rate does not lead to the optimal decay rate of the solution. Motivated 
by using a negative Sobolev space {s > 0) to replace norm, in [lOj we developed 

a general energy method of using a family of scaled energy estimates with minimum derivative 
counts and interpolations among them to prove the optimal decay rate of the dissipative equations 
in the whole space. The method was applied to classical examples such as the heat equation, 
the compressible Navier-Stokes equations and the Boltzmann equation. The main purpose of this 
paper is to applying this energy method to prove the optimal decay rate of the solution to the 
two-species Vlasov-Poisson-Boltzmann system (jl.ip . equivalently, the system (II. 8p for [f,g]. 

Now we define the equivalent instant energy functional by 
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E IIp^/'^h + E ii^^v-^f (1-19) 

h\ + \l3\<N |7|<A^+1 



and the corresponding dissipation rate by 



h\+m<N ^^_20) 

l<|7l<Ar \'y\<N |7|<iV+l 



DECAY OF VPB 



5 



We remark that the definitions of £n and are carefuhy designed to capture the structure of the 
system (11. Sp . Notice that only ||Pi/||^ in <Sjv is not bounded by the dissipation V^- We assume 
that is a sufficiently large integer. Our first main result can be stated as follows. 

Theorem 1.1. There exists an instant energy functional 8]^ [t) such that if £i\f{0) is sufficiently 
small, then the Vlasov-Poisson-Boltzmann system ()1.8p admits a unique global solution [f{t, x, v),g{t, : 
satisfying that for all t > 0, 

ft 



£N{t)+ [ VN{T)dT <£n{0)- 

Jo 

If further, /o G LlH~^ for some s G [0, 3/2), then for all t > 0, 

||A-^/(t)|| <Co, 

and the following decay results hold: 

£N{t)<Co{l + t)-' 

and 



l<k<N 



l<k<N+l 



(1.21) 

(1.22) 

(1.23) 
(1.24) 



Remark 1.2. We remark that by the Poisson equation (jl.SP o. we have 



£. 



N 



|7| + |/3|<Af 



+ ||V,A-ip 



25 



(1.25) 



Remark 1.3. Notice that we do not need to assume that the L^H~^ norm of fo is small. This 
norm is preserved along time evolution and enhances the decay rate of the solution with the factor 
s/2. The constraint s < 3/2 comes from applying Lemma \A.5\ to estimate the nonlinear terms when 
doing the negative Sobolev estimates via . For s > 3/2, the nonlinear estimates would not work. 
We remark that we do not require the L'^H~^ norm of go. 

As we metioned before that the linear homogeneous system of (jl.Sp is decoupled into the Boltz- 
mann equation for / and another system for g and V^^J* looks like the one-species Vlasov-Poisson- 
Boltzmann but with the linearized collision operator £2- However, indifferent from the Boltzmann 
equation, the Vlasov-Boltzmann equation and the one-species Vlasov-Poisson-Boltzmann system, 
this special coupling effect between the Poisson equation and the linear operator £2 gives the dis- 
sipation estimate of not only the microscopic part {I — ^2)9 but also the hydrodynamic part P25 
(it only has one hydrodynamic field, namely, d{t,x)). This is also what we accorded to when we 
define the dissipation Pjy. Then it would suggest that / will decay at the optimal algebraic rate 
as the Boltzmann equation but g (and Va;$) may decay at the exponential rate. These facts can 
be verified easily for the linear homogeneous system of (jl.Sp along our proof of Theorem 11.11 or by 
the Fourier analysis as in [21 [3]. However, for the nonlinear problem (jl.Sp . to control the nonlinear 
terms we need to impose a bit stronger assumption for the initial data. This is our second main 
result which can be stated as follows. 

Theorem 1.4. Under the assumptions of Theorem \l.l\ if additionally £n{^) + ||[/o7 5o]||^ is suffi- 
ciently small, then for all t > 0, 



^ \[v^f{t),V^g{t) 



2<k<N 



+ E ||v''V,$(t) < Co(l + t)-('+^^ (1.26) 

2<k<N+l 
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if additionally <SAr(0) + X^o<fc<7V-i || [^'^/oi ^'^S'o] ||^ is sufficiently small, then for all t > 0, 



0<fc<iV-l 



+ E 

0<fc<iV 



V'=V,$(t) 



< Cqs ^* for some A > 0, 



and 



V^g(t)||^ + ||V^+iV.$(i)|r < Coil + 

E h'fit) '<Co(l + t)-(^+^) /or £ = 3,..., TV- 1, 



i<k<N 



I V'={I - Pi}/(t) < Co(l + /or A: = 0, . . . , iV - 2. 



(1.27) 
(1.28) 

(1.29) 
(1.30) 



We will prove Theorem 11.11 and Theorem 11.41 by the energy method that we recently developed 
in |10) . As there, we may use the linear heat equation to illustrate the main idea of this approach 
in advance. Let u{t) be the solution to the heat equation 

dtu - Au = in R' 

u\t=o = Uq, 

Let —s<£<N. The standard energy identity of (|1.31|) is 

ld_ 
2di 



(1.31) 





2 




2 




+ 


= 




L2 




L2 



Integrating the above in time, we obtain 



V^u{t) 



L2 



< 



V^no 



L2 



(1.32) 
(1.33) 



Note that the energy in (|1.32p is not bounded by the corresponding dissipation. But the crucial 
observation is that by the Sobolev interpolation the dissipation still can give some control on the 
energy: for —s < i < N,hj Lemma lA.31 we interpolate to get 

i+s 

^'+'u{t) 



L2 



V'u{t) 

Combining (fL34l) and (fOSl) (with I 

V'+\{t) 



< \\A-'u{t)\\l+'+' 
= —s), we obtain 



l+l + s 
L2 



L2 



> ||A-"uo| 



1 

1 + 3 



L2 



V^uit) 



1+ 

L2 



e+s 



Plugging (ll.35j) into (jl.32p . we deduce that there exists a constant Co > such that 

1- 



d 
It 



L2 



Co 



1 

' t+s 



< 0. 



Solving this inequality directly, and by (ll.33p . we obtain the following decay result: 

2 / 2_ r< + \ "(^+*) 

£ e+s Ool 
L2 



V^n(t) 



L2 



+ 



+ S 



< Coil + t) 



^{e+s) 



(1.34) 
(1.35) 

(1.36) 
(1.37) 



Hence, we conclude our decay results by the pure energy method. Indifferent from the LP-Lp' 
approach, an important feature here is that the norm of the solution is preserved along time 
evolution and this norm of initial data enhances the decay rate of the solution with the factor s/2. 
Although (|1.37|) can be proved by the Fourier analysis or spectral method, the same strategy in 
our proof can be applied to nonlinear system with two essential points in the proof: (1) closing the 
energy estimates at each £-th level (referring to the order of the spatial derivatives of the solution) ; 
(2) deriving a novel negative Sobolev estimates for nonlinear system which requires s < 3/2 (n/2 
for dimension n). 
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In the rest of this paper, except that we wih cohect in Appendix the analytic tools which will 
be used, we will apply the energy method illustrated above to prove Theorem 1 1 . 1 1 and Theorem 11.41 
in Section [2] and Section [3] respetively. However, we will be not able to close the energy estimates 
at each i-th level as the heat equation. This is caused by the "degenerate" dissipative structure of 
the linear homogenous system of (jl.Sp when using our energy method. More precisely, the linear 
energy identity of the problem reads as 



2dt 



+ 



It is well-known that Ci and £2 are only positive with respect to the microscopic part respectively, 
that is, there exists a cto > such that 



{CiV'f,V'f)>ao V^{I-Pi}/ and {C2V'g,V'g) > ao V^{I - Pa}^ 



(1.39) 



To rediscover the dissipative estimate for the hydrodynamic part [Pi/, P2(7] and the electric field 
Va;$, we will use the linearized equation of (jl.Sp via constructing the interactive energy functional 
Gf between and V^^^ of the solution to deduce 



dt 



+ 



V^+ipi/ 



L2 



+ 



V'P2<7 



L2 



+ 



L2 



+ 



< 



V^{I-Pi}/,V^{I-P2}5 



L2 



+ 



v^+Hi-Pi}/,v^+Hi-P2}<? 



2 

L2 ' 



2 

L2 



£+1 J 



(1.40) 



This implies that to get the dissipative estimate at each level for the missing part in the standard 
energy inequality it requires us to do the energy estimates (jl.SSp at two levels (referring to the order 
of the spatial derivatives of the solution). To get around this obstacle, the idea is to construct some 
energy functionals £i{t)^ < H. < N — \ (less than — 1 is restricted by ()1.40p ). 

E ||^'/(^) 

i<k<N 



L2 



which has a minimum derivative count of and we will derive some Lyapunov-type inequalities for 
these energy functionals in which the corresponding dissipation (denoted by T)i{t)) can be related 
to the energy similarly as (jl.SSp by the Sobolev interpolation. This can be easily established for 
the linear homogeneous problem along our analysis, however, for the nonlinear problem (jl.Sp . it is 
much more complicated due to the nonlinear estimates. This is the second point of this paper as in 
[lOj that we will use extensively the Sobolev interpolation of the the Gagliardo-Nirenberg inequality 
between high-order and low-order spatial derivatives to expect to bound the nonlinear terms by 
^J£N{t)V|>{t) that can be absorbed. But this can not be achieved well at this moment and we will 
be left with two extra terms: one term is related to a sum of norms of V^^ and the other term is 
related to a sum of velocity- weighted norms of [/,(?], as stated in (j2.78p . Note that when taking 
£ = 0, 1 in (|2.78p . we can absorb these two unpleasant terms and then we get Theorem 1 1 . 1 1 after we 
complete the negative Sobolev estimates. While for ^ > 2, we need to assume the weighted norm 
of the initial data. With the help of these weighted norms, we will succeed in removing the sum of 
velocity- weighted norms of [/, 5] from the right hand side of (I2.78P to get (I3.29P in which we can 
take 1 = 2. On the other hand, we can revisit the equations (|1.8P o- (|1.8P n to deduce a further energy 
estimates (I3.2ip for g and V^^ which kills the sum of norms of V^*^ in the right hand side of (I3.29P 
to get (|3.33p in which we can take £ = 3,...,A^ — 1. This energy estimates (|3.2ip also implies the 
exponential decay of g and Va;$. Finally, to estimate the negative Sobolev norm in Lemma 12.41 
we need to restrict that s < 3/2 when estimating A~'^ acting on the nonlinear terms by using the 
Hardy-Littlewood-Sobolev inequality and also we need to separate the cases that s G [0, 1/2] and 
s G (1/2, 3/2). We remark that it is also important that we use the Minkowski's integral inequality 
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to interchange the order of integrations in v and x in order to estimate the nonhnear terms and 
that we extensively use the spHtting / = Pi/ + {I — Pi}/ and g = V^g + {I — ^2}g- Once these 
estimates are obtained, Theorem 11.11 and Theorem 11.41 fohow by the interpolation between negative 
and positive Sobolev norms similarly as the heat equation case. 

To end this introduction, we want to emphasize that the results of the two-species Vlasov- 
Poisson-Boltzmann system in Theorem 11.11 and 11.41 reveal the essential difference when compared 
to the one-species Vlasov-Poisson-Boltzmann system or the two-species Vlasov-Boltzmann system. 
For the one-species Vlasov-Poisson-Boltzmann system |2], the electric field decays at the optimal 
algebraic rate; For the two-species Vlasov-Boltzmann system [22j the disparity between two species 
of particles decays at the optimal algebraic rate as the Boltzmann equation. Besides this, our 
achievement heavily relies on the special cancelation property between two species that our refor- 
mulation (ll.Sp displays well which gives the dissipative estimates of the norm of electric field. 
This cancelation phenomenon was also observed in [l_9j for the study of the two-species Vlasov- 
Maxwell-Boltzmann system. The natural generalization of this paper is to using our energy method 
to revisit the decay rate of the two-species Vlasov-Maxwell-Boltzmann system [23] . Another inter- 
esting application is to revisit the one-species Vlasov-Poisson-Boltzmann system system in which 
the norm of electric field is not included in the dissipation, and we expect to investigate it in 
the future study. 



2. Energy estimates and proof of Theorem 11.11 

2.1. Basic energy estimates. In this subsection, we will derive the a priori nonlinear energy 
estimates for the system (jl.Sp . We shall first establish the energy estimates on the pure spatial 
derivatives of solutions. First of all, notice that the dissipation estimate in ()1.39p is degenerate, 
and it only controls the microscopic part [{I — Pi}/, {I — P2}5]- Hence, to get the full dissipation 
estimate we shall estimate the hydrodynamic part [Pi/, P25'] and the electric field V in terms 
of the microscopic part, up to some (small) error terms. This will be done in the following lemma. 

Lemma 2.1. If £j\[{t) < 6, then for k = 0, . . . , N — 1, we have 



dG) 
dt 



+ C 



< 



v^i-Pi}/ + V^+Hi-Pi}/ 

2 



+ 6' 



+ 



(2.1) 



and 



d_ 
It 



< 



+ 



V'=+iP25 



+ 



VHl-P2}<7 



+ 



+ 5' 



(2.2) 



Here Gh(t) and G^{t) are defined by (|2.18p and (j2.23p respectively which satisfy the estimates 



G) 



< 



+ 



and 



g: 



< 



V^g 



+ 



+ 



(2.3) 



Proof. As in |51 HH 1^ [2t [5], the proof is based on the careful analysis of the local conservation laws 
and the macroscopic equations which are derived from the so called macro-micro decomposition. 
Our contribution is to carefully estimate the nonlinear terms so that we may close the energy 
estimates at each ^-level in our weaker sense. 
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First, multiplying (jl.Sp ^ by the collision invariants l,v, \v\'^/2, (jl.SP n by 1, and then integrating 
in t; G M^, we get the local conservation laws 



dt Fdv + V^- vFdv = 0, 

dt / vFdv + Vj,- / v(^vFdv -Vr,<^ I Gdv = 0, 



dt [ hv\'^Fdv + V^- [ hv\'^vFdv -V^^ ■ [ vGdv = 0, 



dt Gdv + V^- vGdv = 



Plugging F = // + ^ (Pi/ + {I - Pi}/) and G = ^{^29 + {I - ^2}9) into ([231), and using the 
representations (I1.13P and (I1.15p . we obtain 

(9ta + • 6 = 0, 

+ V, • (a + 4c) + V, • ^({I - Pi}/) = d, 
dt{2>a + 12c) + V, • (56) + V, • B{{1 - Pi}/) = V,«> • V{{1 - Pa}^), ^^'^^ 
atd + V,-P({I-P2}5) = 0, 

where we have defined the moment functions A = (^jjOsxSj ^ = (■Si, ;S2, iSa) and V = (2?i,P2)^3) 
by 

Aij{h) = {h,ViVj^/JT), Bi{h) = {h,\v\'^Vi^/Ji) and T>i{h) = {h,Vi^/Jl). (2.6) 
Notice that (|2.5p ^ and (j2.5P o implies 

5tc + Jv, • 6 + i-V, • B{{I - Pi}/) = -^V,«> • V{{I - FM. (2.7) 
12 12 

Also, the Poisson equation (jl.SP n reads as 

= d. (2.8) 
Next, plugging / = Pi/ + {I - Pi}/ and 5 = P25' + {I - ^2)9 into ([TH]), we obtain 
dtPif + ^ • V^Pi/ = - Pi}/ + £1 + OTi, 

5tP25 + • VxP29 - • = -dt{I - P2}9 + ^2 + ^2, 
where the linear terms £1 and £2 are denoted by 

£i = -^^-V,{I-Pi}/-£i{I-Pi}/, £2 = • V,{I - P2}5 - £211 - P2}5. (2.10) 

Motivated by [21 [3], we will not use all of the macroscopic equations but only use the evolution of 
some moments of the microscopic part that appeared in (|2.5p . To simplify the notations, we define 
two more moment functions A = (Aij) 3x3 and B = {Bi,B2,Bs) by 

Aij{h) = A^J{h) - {h, B,{h) = ^ {Bi{h) - 5Vi{h)) . (2.11) 

Applying Aij{-),Bi{-) to i\2.9^ ^ and I^i(-) to (l2.9P n respectively, we get 

dt[Aj{{i - Pi}/) + 4c5,j] + d,bj + djbi = A,(£i + oil), 

dtB,{{I - Pi}/) + 2dic = B,{£i + ail), (2-12) 
dtV,i{I - P2}5) + did - di^ = A(£2 + ^2). 



(2.9) 
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Notice that for fixed i, one can deduce from (|2.12p ^ that 



dt 



^ a, - Pi}/) + ^d,Au{{i - Pi}/) 

j 

\ Yl + Oil) - ^ 9, + 



A^^6,; - diibi 



(2.13) 



Bypassing the nonhnear couphng terms, the local conservation laws ()2.5P ]^- p.5P o. ()2.7p and the 
macroscopic equations p.l2P ]^- ()2.12P n. ()2.13p on the macroscopic coefficients a,b,c are decoupled 
from those on d and V^<I>, and they have the same structure as the pure Boltzmann case. So by 
these equations and following the proof in [8] or [U O O [H] , we can deduce the dissipative estimates 
of o, b, c: for I7I = k with A; = 0, . . . , iV — 1, 



^ |^a'^^9,A,({i-Pi}/) + ^a^a,A.({i-Pi}/),5^6,j +i 



< e ( Wd^V^af + \\d''V,.cf + \\d^ (V,$d)f 



+ c, \\d^{i - Pi}/f + ||v,5^{i - Pi}/f + \\dm,. 



(2.14) 



dt 



d^Bii{i-Pi}f),dWic) + ||a^v,.cf 

< e (||a^V,6f + \\d^' {V,^ ■ V{{I - P2}5))ir 



(2.15) 



+ Ce - Pi}/f + ||v,a^{i - Pi}/f + I 

^(a^6,9^V,a) + i||a^V,af 

< c (||a^v,.6f + iia^v.cf + ||v,a^{i - Pi}/f + Wiv,^ , 



(2.16) 



Here is defined as (9Ti,C) with C is some linear combination of [.^,Vi^/J2,ViVj^/Jl,Vi\v\'^ ^/Jl], 
etc. We multiply ()2.14p and (j2.15p by a constant M > and then sum up them as well as (j2.16p . 
We first fix M to be sufficiently large such that the first two terms in the right hand side of (j2.16p 
can be absorbed, and then for such fixed M we further let e > sufficiently small such that the 
first term in the right hand side of (|2.15p and the first two terms in the right hand side of (j2.14p 
can be absorbed. Hence, we obtain 



'^-G'lit) + C 



V'^+^Pi/ 



dt f 

where G^{t) is defined by 



< 



+ 



V^U-Pi}/ 
V^(V,$d) 



+ 



+ 



v'^+Hi-Pi}/ 



+ 



(V,$-P({I-P2}5)) 



(2.17) 



G}{t):= Y, \ M(a^j;a,^,({i-Pi}/) + i9^a,A,({i-Pi}/),a^6,) 
i7i=n V i / 

+M (a^s.({i - Pi}f),d''dicj + {dn, d^v^a)} . 



(2.18) 



We now focus on the derivation of the dissipation estimates of d and V^;^. However, we shall 
estimate in the same spirit with a bit more attention on the electric field. Applying d'^ with I7I = k 
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to (j2.12P 3. multiplying the resulting equations by d'^{did — and then integrating by parts over 
G M'^, by the Poisson equation (j2.8p . we get 



= 2{d^di^, d^did) - - P2}5), d^{did - 

+ {d^'Vi{2.2 + 'yi2),d\did-di^)) 



-2\\d"<d\[ 



d_ 

di 



id^Vii{I-P2}9),d^d,d-di'^)) 



+ {d^Vi{22 + 'yi2),d^did-di<^>)) 
< -2 wd^df - ^ - P2}9),d'y{d,d - dm 

+ £ (\\dWtdf + lia^StV^cDf + \\d^V,df + ||9^V,.^'| ^ 

+ Ce - P2}5f + ||V,.9^{I - Psifff + \\dmo\\ f 



(2.19) 



It is crucial that it follows from (j2.5p ^ and (j2.8p that the follows hold: 

Wd^dtdf <c\\v.d^{i-P2}gf 



(2.20) 



and 



(2.21) 



Hence, by (|2.20p - ()2.2ip and the Poisson equation (j2.8p again, we deduce from ()2.19p that, by taking 
e > sufficiently small. 



d_ 

di 



V'^P^g 



+ 



v'+'P2g 



+ 



< 



V^II-Pajg + V'^+Hl-P2}5 



2 



+ 



where Gg(t) is defined by 



Ggit) ■■= E (5"A({I-P2}5),5^(5.d-a,$)). 

i7|=A: 



(2.22) 



(2.23) 



It is clear that G^{t) and Gg{t) satisfy the estimates (j2.3p . and it remains to estimate the 
nonlinear terms in the right hand side of (j2.17p and (j2.22p . The main idea is that we will care- 
fully adjust the index of spatial derivatives when estimating the nonlinear terms so that they can 
be bounded by the right hand side of (j2.ip or (j2.2p . We begin with the estimate of the term 
||V''ai2,||||^ := E|7|=fc \\{d^^2X)f- First, by the estimate ([02]) of Lemma|X8]and the fact that 
^ decays exponentially in u , we have 



J] IKa^r(5,/),C)f < E E ||(r(5^^5,a^-^V),C)| 

|7|=fc |7|=A:7i<7 

l7i|<fc 



(2.24) 
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By Holder's inequality, Minkowski's integral inequality (|A.14p of Lemma lA. 61 the Sobolev interpo- 
lation of Lemma IA.2I and Young's inequality, we obtain 



< 






^7^^-1711 f 


< 



















il- 



< 5 



l7ll + l 



+ 



l7ll + l 
fc + 1 



l7ll + l 



l7ll + l 
fc + 1 



Here a comes from the adjustment of the index and is defined by 

1 fc-|7i| _ /I a\ |7i| + l (I k + l 

3 3 \2 3yA; + l 

k + l k + l ]\ 

=^ a 



X 1 



l7il + l 
k + l 



k+l k+l N 
2(|7i| + l) -^-y 



Hence, we have 



^ Hiding, f),Of< 6' ( v^+v + 
Similarly, we may apply the same arguments to obtain, by Lemma lA.ll and Lemma lA.21 

\'Y\=k |7|=fc7l<7 



+ 



<6' 

and with the additional integration by parts over t;-variable to have 











2 




+ 



^ ll(a^(v.ci>-v./),c>f = IK5^(v.<i>/),v.c>f 

|7|=fc |7|=fc 



|7|=fc 7i<7 
ir2 / llvrfe+l/ 



<,5^ V«+^(/,V,C> 



+ 



+ 



Thus, summing up the estimates (|2.27p - (|2.29p . we have 



< 5' 



+ 



Then plugging the estimate (I2.30p into (I2.22p . since 5 is small, we obtain (12. 2p . 
Similarly, we have 



< 6' 



+ 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



(2.29) 



(2.30) 



(2.31) 
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However, we will also use the same argument above to estimate the remaining two terms in the 
right hand side of (j2.17p to obtain 

2^ 



l7l=fc 







2 
+ 




<5'( 




2 

+ 





and 



^ \\d^ (v,.<i> • v{{i - P2}g))f < ( v^+ip({i - P2}9) 

h\=k 



+ 



<6' 



+ 



Then plugging the estimates (I2.3ip - (j2.33p into (I2.17p . since 5 is small, we obtain (12. ip . 

Now we consider the energy estimates for the pure spatial derivatives of the solution. 
Lemma 2.2. If £iy{t) < 6, then for k = 0, . . . , N — I, we have 
d 



dt 



+ 



+ c 



<6 



V'^{I-Pi}/,V'''{I-P2}<7 

2 



+ 



2<e<N 

and for k = N , we have 



+ ||[v^{i-Pi}/,v^{i-P2}5]|[) ; 



d_ 
dt 



[V^/, V^5] II' + II V^V.c^ll^j + C II [V^{I - Pi}/, V^{I - P2}g] li: 
<<5(||[v^/,v^5]||'+ E ||vV.ci,|p+ Y ||[v^{I-Pi}/,V^{I-P2}5 

y 2<i<N i<e<N 

Proof. The standard with I7I = A; energy estimates on (jl.Sp ^ and (II.8P 2 give rise to 



and 



Ijt + ^^^^ = (dmud^) := /i + I2 + /3 



Ijt " ^^^^^'■^ ■ + (^29^5,9^5) = (9^^2,9^5) := Ji + J2 + J3. 



We first estimate the left hand side of ([236]) -([237]). By (l23])4 and ([23]), we have 
-(a^V,«> • v^, d^g) = V, • V{{1 - P2}5)) 

The estimate (|A.19p of Lemma IA.7I implies 



and 



iC,d^f,d^'f)>ao\\d^{I-Pi}f\\l 



iC2d^g,d^g)>ao\\d^I-P2}g\\: 



(2.32) 



(2.33) 



□ 



(2.34) 



(2.35) 



(2.36) 
(2.37) 

(2.38) 

(2.39) 
(2.40) 
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(2.41) 



Next, we turn to estimate the right hand side of (j2.36p - (|2.37p . First, by the cohision invariant 
property and the estimate ()A.23P (with r] = 1/2) of Lemma lA.81 we obtain 

Ji := {d^T{gJ),d'^g) = ^ C7^^ {r{d''^ g,d''-^^ f),d^I - I'M 

< \\^-'/Md'''9,d''-^'f)\\\\u'/^d^{I-F2}g\ 

< E |||Vl^^l5|2|V^'-l^il/|, + |Vl^^l5|.|V'=-l^^l/|2||||v^'{I-P2}5 

71 <7 

For the first term in the right hand side of (j2.4ip . we use the splitting / = Pi/ + {I — Pi}/ to 
have 

||Vl^^l5|2|V'=-lT'^l/U|| < |||VlT'^l5|2|V'=-lT'^l/|2|| + |||Vl^^l5|2|V^-l^il{I^ 
:= Jii + Ji2- 

For the term Jn, if A: = 0, . . . , — 1, it has been already bounded in ()2.25p as 



(2.42) 



Jn <6 



+ 



(2.43) 



while for k = N, we then take the supremum in x on the factor with the less number of derivatives 
to have, for instance if |7i| < by Lemma lA.61 and Lemma lA.21 



Jn < 



VlTll, 



<l|v".r-^||v^||^||/||¥||v-/||-^ 



(2.44) 



<5(||V^/|| + ||V^5||), 
where we have denoted a by 



3 



1 a 

2 ~ 3 

3A 



X 1 



M 

A^ 



+ 



1 N\ |7i| 



2 3 



N 



(2.45) 



a 



< 3 since |7i| < A/2. 



2(A^-|7i|) 

For the term J12, note that we can only bound the z^- weighted factor by the dissipation, so we can 
not pursue as before to adjust the index. Noticing that {I — Pi}/ is always part of the dissipation. 
If ItiI < k — 2 (if A; — 2 < 0, then it's nothing in this case, etc.) and hence A; — |7i| > 2, then we 
bound 



J12 < vi^ii 



9\ 



V^-|tiI{i_p^}/ 



<5 E ||v^{I-Pi}/ 

2<1<N 

and if I71I > A; — 1 and hence A: — |7i| < 1, then by Sobolev's inequality, we bound 



(2.46) 



J12 < ||vi^^i<7|| IllV^-i^iiji-Pi}/!, 

<5 E ||v^{I-Pi}/"' 

1<£<3 



< 


V\'<^\g 











,l/2yfc-|7i|{l_p^}J 



Hence, we have 



J^^^^ E ||vni-Pi}/ 



(2.47) 



(2.48) 



l<^<Af 
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Thus, we complete the estimate of the first term in the right hand side of ()2.4ip . Applying the 
same argument (exchange / and ^f) to the other term, we conclude that for A; = 0, . . . , — 1, 



+ 



v 



HI-P2}5|[+ Yl ||[v^{I-Pi}/,V^{I-P2}5]||j 

l<i<N "J 



and for k = N, 



Ji<<5(||[V^/,V^5]||+ E ||[v^{i-Pi}/,v^{i-P2}5]f I . 

y i<i<N "J 

Similarly, we have that for A: = 0,...,A — 1, 

h:={d'rr{f,f),dy)<sl V'^+V + V^{I-Pi}/ + ||v^{I-Pi}/ 



(2.49) 
(2.50) 



i<e<N 



and for k = N, 



/i<W||V^/||+ Y ||v'{i-Pi}/ 

\ l<^<Af 

Next, for the second term J2 in the right hand side of (|2.37p . we first have 

J2 := (5^(V,cI> • vf),d^g) = ^ C:;Hd^^V,^ ■ vd''-^^ f g) 

71 <7 

< |||vi^iiv,<i>||v'=-i^ii/|,||||v* 

71 <7 

Applying the same arguments that for (|2.42|) . we obtain that for A; = 0, . . . , A — 1, 



(2.51) 



(2.52) 



(2.53) 



+ 



+ 



+ Y ||v^{i-Pi}/ 



l<t<N 



(2.54) 



(2.55) 



J2<5 
and for k = N, 

^2< J||V^/|| + ||V^V,<i>|| + ||V^5|L+ Y ||v'{i-Pi}/ 

\ l<^<Af 

For the second term I2 in the right hand side of (|2.36|) . we again use the splitting to have 
h := (5^(V,<I> • vg),dy) = Y qH5^^V,<i> • vd'y~^'g,d'yf) 

71 <7 

= Y C:!; {id^'V^^-vd'^-^''g,d'^Pif) + {d'^'V^^-vd^-'^'g,d^{I-Pi}f)} (2-56) 

71 <7 

:= hi + -^22- 

Notice that the term I22 can be bounded as that for ()2.42p and we have that for A; = 0, . . . , A — 1, 



I22 < 6 



+ 



+ 



V^I-Pi}/ + V ||V^{I-P2}5 



(2.57) 



l<^<Af 
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and for k = N, 

^22<5||||V^<7|| + ||V^V,f|| + ||V^{I-Pi}/||^+ Yl ||V'{I-P2}5 
While for the first term I21, we have 



i<^<Af 



l2i< I |Vl^ilV,<i>| V^-I^^l<7 iV^bdx 



(2.58) 



(2.59) 



For < k < N — 1, hy Holder's inequality, Minkowski's integral inequality ()A.14p of Lemma lA.61 
and Sobolev's inequality, we have 



/21< 

< 



LI 



< 



LI 



If 7i = 7, then in this case we have 



LI 



< 6 



+ 



+ 



(2.60) 



(2.61) 



If I71I <k — \ (it is nothing if A; — 1 < 0), then in this case by Lemma [A . 1 1 and Lemma [A. 21 we have 



~ II 11^ 

WnW k 



l7ll 



< 5 



+ 



V^g 



+ 



where we have denoted a by 
3 3 



l7i| 



(2.62) 



k 



2 3 / k 



"= 2(fc-|7i|) -^ since |7i|<fc-L 



(2.63) 



Hence for 0<A;<A^ — 1, we obtain 



+ 



+ 



+ 



(2.64) 



Now for k = N, we may use the arguments as in ()2.44p to have 

/2i<<)-(||vM|%||v^v.<i>||%||v^/||^ 

Hence, in light of (l2371) - ([238]) and we may conclude that for k = 0, 



(2.65) 



,A^- 1, 



+ 



V^'{I-Pi}/ + Yl ||V^{I-P2}5 



+ 



V^g 



(2.66) 



1<^<JV 



and for k = N, 



/2<5 I ||[V^/,V^5]|r + ||V^V.cI>|| + ||V^{I-Pi}/||^+ ||V^{I-P2}5 



i<e<N 



(2.67) 
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Now we turn to the third term J3 in the right hand side of (j2.37p . There is one worst case 
involving {k + l)-th derivative. We should estimate this {k + l)-th order derivative term together 
with the similar term stemming from the third term I3 in (|2.36p to be canceled by the integration 
by parts over T;-variable. Hence, we estimate them together to have 

/3 + J3 := (a^(v.<i> . v.f),d^g) + (a^(v.<i> • v.g),dy) 

+(a^w,<i> • v,a^-^i5, 5^Pi/) + (5^^ v,«> • v,a^-^i5, - Pi}/)} • 

(2.68) 

For the third term in the right hand side of (|2.68|) . we integrate by parts in v and use the fast that 
Pi/ decay exponentially in v to have 

(a^^v^cD • v„a^-^i<7,5^Pi/) = -(a^^v,«>a^-^ig, v.a^Pi/) 

(2-69) 



iVVbdx. 



This term appeared in (|2.59|) that has already been bounded. Similarly for the first term, 



< 



V 



fc-|7i 



/ 



\V''g\2dx, 



(2.70) 



This term can be bounded similarly as the first term in (|2.42p . While for the remaining two terms 
in (j2.68p . note that we can only bound the factors involving the velocity derivative by the energy 
(not dissipation; otherwise, our method would not work.), so we can not pursue as before to adjust 
the index. Thus, we may bound, for instance the third term, if 1 < |7i| < N/2, 



V.V'^-lTil/ V'={I-P2}5 



<5 



\2<^<[f ]+2 



+ 



and if 1 71 1 > N/2 



/ 



v^v'^-i^ii/ 



^^{I-Ps}^ 



V'={I-P2}5 



(2.71) 



< 6 



Y ||W,a> + V'={I-P2}ff 

\iq-]<i<N 



(2.72) 



Hence, we may conclude that for < A; < — 1, 



^3 + J3<M l|V'+V|| + 



+ 



V''V,,.^> 



+ 



+ Yl 11^'^-^ + [vHl-Pi}/,V'={I-P2}5 



2<e<N 



(2.73) 
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and for k = N, 



<5 ||[V^/,V^5]f+ E IK'^--^ +||[V^{I-Pi}/,V^{I-P2}<7] 



2<1<N 



(2.74) 



Consequently, plugging the estimates (|2.38p - (|2.40p into the left hand side of (|2.36p - (|2.37p and 
then bounding the right hand side by (ITiQil . (l23ni . (l23il) . (I2:66D and ([27731) for < /c < iV - 1, 
and bounding it by (12301) . (l232|) . (1135]) . (I237I) and ^Ol^ for /c = TV, since 5 is small, we obtain 
(I23ill and ([235]) respectively. □ 



Until to now, we may conclude our energy estimates of the pure spatial derivatives as follows. 
Let < a < N — 1 and assume the a priori estimates that Sis[{t) < 6 for sufficiently small 5 > 0. 
Then summing up the estimates (I2.34p of Lemma 12.21 from A: = ^ to — 1 and adding the resulting 
estimate with the estimate ([2.35p . by changing the index, we obtain 



d \ - 

dt ^ 

e<k<N 



E ||[v'={I-Pi}/,V'={I-P2}5 

^ £<k<N 

E ||v^-P25f + E 11^'^-^ 

fc<Af-l e<k<N 

+ E ||v^V,<i>'+ E ||[v'{I-Pi}/,V^{I-P2}5]|[) 



<C26\ E II [vV,^^^ 

,e+l<k<N 



(2.75) 



2<k<N 



On the other hand, we sum up the estimates (I2.ip - (l2.2p of Lemma l2.ll from k = I to N — 1, hy 
changing the index, to obtain, since 6 is small. 



e<k<N-l 





V^Pi/ 


'+ E 


\i+l<k<N 




e<k<N 



+ E 

e<k<N+l 



<C, Yl ||[v^{I-Pi}/,V^'{I-P2}5 



i<k<N 



(2.76) 
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Then, multiplying (j2.76p by a small number /? > and then adding the resulting inequality with 
(j2.75p . we obtain 



d_ 

dt 



i E (II [vv,^^,] 

\e<k<N ^ 



+ 



^ e<k<N-l J 



+ (Ci-C4/3) ||[v'={I-Pi}/,V'={I-P2}<7]|| 

£<k<N ^ 

+ ^3/3 1 llv^Pi/llV IK'P29|f+ Y ll^'^v-'^ 

y+l<fc<Af l<k<N l<k<N+l 

<c25l Y Wl"^' fy'a] '+ Y IK'Pa^V Y 11^'^-^ 



,e+l<k<N 



(2.77) 



Y ||v'=v,ci>||V Y ||[v'={i-Pi}/,v*^{i-P2}5 



2<k<N 



l<k<N 



We define £i{t) to be the expression under the time derivative in (j2.77p . We may now fix (3 to 
be sufficiently small so that (Ci — 2Cif3) > and that £i{t) is equivalent to 



E II [v/, 



V'9 



i<k<N 



+ Y ||v*^v^^ 

^<A;<Af+l 



where we have used the estimates (j2.3p and the Poisson estimate. On the other hand, since (3 is 
fixed and 6 is small, we can then absorb the first three terms in the right hand side of (j2.77p to 
deduce that for i = 0, . . . , N — 1, by adjusting the constant in the definition of £i{t)^ 



j^£i+ Y ll[vni-Pi}/,v'={i-P2}<?]||^ 

i<k<N 

+ Y iiv'Pi/f+ Y ||v'P2<? V Y 



l+l<k<N 



e<k<N 



e<k<N+l 



(2.78) 







'+ Y [V''{I-Pi}/,V^'{I-Pi}<7] 


\2<fe<Af+l 




l<k<N 



Now we turn to the energy estimates on the spatial-velocity mixed derivatives of the solution. 
First notice that for the hydrodynamic part [Pi/, P25'], 



djPif,d],P2g <||[a^Pi/,a^P25]|| 



(2.79) 



which has been estimated in (j2.78p . it suffices to estimate the remaining microscopic part 



a^{i-Pi}/,a^{i-P2}5 



for I7I + |/3| < iV with |/3| > 1 (and hence I7I < iV - 1) 
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Lemma 2.3. If £^{1) < 5, then we have 
Jt E \\\dl{l-^i}f.d}{l-V^}g\ \C \\\d}{l-^i}f.d}{l-V^}g 



|7| + |/3|<]V 
I/3|>1 



|7| + |/3|<JV 
I/9|>1 



< 



E II [v/. 



V'9 



l<^<Af 
+ 



E 



0<^<Af-l 



(2.80) 



||[v^{I-Pi}/,V^{I-P2}5 



(2.81) 



Proo/. We take dj (with |/3| = m > 1) of the equations ([LiDi-dLSia to get 
dtdlil - Pi}/ + V ■ V,a^{I - Pi}/ + d}C^{l - Pi}/ 
+ dtd^-p^f + • V,a;^Pi/ + C^-^dp^v ■ VccdJ_^J 

= dj (r{f, /) + ^V^<^ ■ vg - ■ , 
dtdj{l - P2}9 + V ■ V^dJ{I - P2}9 + dJC2{I - P2}9 

+ dtdJP^g + V • V^d^P^g + C^p'd^.v ■ V^dl_p^g - d^V^^ ■ d^(v^) ^^.82) 



= dj ^r{g, f) + -V,^, . vf - • V,/ j . 

We ihustrate only the estimate on {I — P2}g and the other one, {I — Pi}/, can be estimated in 
the same way. Taking the inner product of (I2.82P with dj{l — P2}g, we obtain 

i| ||a^{I - P2}gf + {d}C2{l - P2}g, dj{l - P2}g) 

+ {dtd}P2g + V ■ V^dJP2g + C^'dp.v ■ V.d1_f,^g - 9^V,<I> • dp{v^), dj{l - Pa}*?) (2.83) 

= (^dj (rig, /) + iv.c^ • vf - V,^ ■ V,/) , dj{l - P2}g^ . 
By the hnear estimate ()A.2ip in Lemma \A.7\ we know 

[dJC2{I - P2}g, dj{l - P2}g) > \ \d}{l - P2}5|[ - C liani - P2}<7||^ • (2.84) 

We now estimate the second line in (j2.83p . For the first two terms, by the local conservation 
laws (j2.5P | and Cauchy's inequality, we get 



(2.85) 



dtdlP2g + v ■V^,d}P2g,dl{l-P2]g] < {\\dtd^ d\\ + \\V ^,d"< d\\) dl{l-P2]g 



1 

< — 
- 12 



dl{i-P2]g +c||v,a^g||^ 



For the third term, we use the sphtting to have, since |/3 — /3i| = m — 1, 
yH-V^d}_p^g,d}{l-P2}g] 



< 



1 

< — 
- 12 



[d^^v ■ V^d;_^P2g, d}{I - P2}g) I + I [d^'v ■ V,a^_^JI - P2}g, dj{l - P2}g 
d]^{l-P2}g\c\\V,d^'gf + C Y ||5^{I-P2}5 



(2.86) 



l7l + l/3|<iV 
|/3|=m-l 
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d'^v,'^-dp{v^),dm-P2}g 



< 



1 

12 



- P2}5 



(2.87) 



Now we turn to the third hne in (j2.83p . First, by Lemma 7 of [19] (clearly, without taking the 
time derivatives), we have 

(^d}T{gJ),d}{I-P2}9) < VS^V^. (2.88) 
Then it remains to estimate the last two terms related to the electric field. We do the splitting 

(^dj Qv.cD . vf - V,<l> ■ V,/) , dj{l - P2}g^ 

= l^dj Qv,.<i> • vPif - V,.«> • V.Pi/^ , d^il - P2}9^ (2.89) 

+ (dj Qv.cD . y{i - Pi}/ _ v,^ ■ V,{I - Pi}/) , dj{l - Ps}^) . 

Since the hydrodynamic part is not affected by the velocity derivative and the v factor as noted in 
(j2.79p . hence the first term in (|2.89p can be bounded by \J£nT^n- The argument is similar as that 
in Lemma [2. 2 1 but is a bit simpler since we do not adjust the index. The main concern is that there 
is one worst case involving (A^ + l)-th derivative in the second term coming from the V^^-Vlasov 
term. We should estimate this {N + l)-th order derivative term together with the similar term 
stemming from the equation (|2.8ip to be canceled by the integration by parts: 

dZ (V,<i> • V4I - Pi}/) , dl{l - P2}g] + (dl (V,<i> • V^I - Pi}5) , dm - P2}/) = 0. (2.90) 



After this cancelation, we then can bound the second term in (j2.89p by \/£nT^n- Plugging these 
estimates and (|2.84p ~ (j2.88p into (j2.83p . and doing the same estimates for (|2.8ip with respect to 
{I — Pi}/, summing up I7I + < N with |/?| = m > 1, we obtain 



1 d 
2di 



\\\^J{I-P^}f,^],{I-P2}9\ +7 E ||k{I-Pi}/,5^{I-P2}5 



7| + |/3|<JV 
|/3|=m 



< 



T. II [v/, 



KKN-1 



+ 



E 



0<£<N-1 



\l\ + \P\<N 
/9|=m 

2 



+ E ||[v^{I-Pi}/,V^{I-P2}5 



0<^<Af 



+ E \\\dJ{I-Pi}f,dJ{l-P2}9 



|7| + |/3|<JV 
|/31 = m-l 



A simple recursive argument on (j2.9ip with the value of m gives (j2.80p . 



(2.91) 

□ 



2.2. Negative Sobolev estimates. In this subsection, we will derive the evolution of the negative 
Sobolev norms of the solution. In order to estimate the nonlinear terms, we need to restrict ourselves 
to s G (0,3/2). We will establish the following lemma. 



Lemma 2.4. If £]\f{t) < 6, then for s G (0, 1/2], we have 



I ||A-V||' + C \\A-'{I - P^}f\\l < (||A-7|| + 1) Vr,; 



and for s £ (1/2,3/2), we have 
d 



^ ||A-/||^ + C ||A-^{I - P,}f\\l < (II A-^/|| + 1)Vn + \\ft+' IIV/II 



l3-2s 



(2.92) 
(2.93) 
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Proof. Applying A ^ to (jl.SP ;^ , and then taking the inner product with A we have 
^ ^ lA-7||' + ^o||A-^{I-Pi}/||' 



2 dt 



(2.94) 



We win estimate the right hand side of ()2.94p term by term. For the first term, by the colhsion 
invariant property, we have 



(A-^r(/, /), A-v) = (A-^r(/, /), A-^{i - p}/) 



< 



A-Uu~-2Tif,f)) z.2A-^{I-P}/ 



(2.95) 



< C 



A-(^-^r(/,/)) +^||A-^{i-P}/||;. 



To estimate the right hand side of ()2.95p . since < s < 3/2, we let 1 < p < 2 to be with 
1/2 + s/3 = 1/p. By the estimate ()A.12p of Riesz potential in Lemma \A.5\ Minkowski's integral 
inequality ()A.14p of Lemma lA. 61 and the estimate ()A.23P (with rj = 1/2) of Lemma \A.8\ together 
with Holder's inequality and the splitting / = Pi/ + {I — Pi}/, we obtain 

2 



A-^ [i^-^TifJ) 

^'5r(/,/) 



< 



< 



^"^r(/,/) 



< 



< 



2I/U1IL? 



(2.96) 



< 



[I-P}/t+ll/ll)- 



We bound the first term in ()2.96p as, since 3/s > 2, by Sobolev's inequality, 

[I - P}/L < \\f\\L?.m - ^ ^ - P}/IL 



3 

LyLx 



(2.97) 



While for the other term in (j2.96p . we shall separate the estimates according to the value of s. If 
< s < 1/2, then 3/s > 6, we use the Sobolev interpolation and Young's inequality to have 



3 11/11 <||V/f+^/2 11^2^ 



2 .||l-s/2 



<5(l|v/|| + ||vV||); 



(2.98) 



and if s G (1/2,3/2), then 2 < 3/s < 6, we use the (different) Sobolev interpolation and Holder's 
inequality to have 



<ii/ir"'/'iiv/f/'- 



(2.99) 



For the second term in (|2.94p . we do the splitting / = Pi/ + {I — Pi}/ to have, similarly as in 

(A-^ {V,^.vg),A-'f) 

= (A-^ (V^cD . yg) , A-^Pi/) + (A-^ (V,.<I> • vg) , A-^{I - Pi}/) 

< ||A-^ {\VM9\2)\\ ||A-7|| + ||A-^ (|V,<I>||5|.)|| ||A-^{I - Pi}/||^ 

< \\\vM9\2\\li ||a-^/|| + \\\vM9\u\\l^^ ||a-^{i - Pi}/L (^-^^^^ 

<||V,<I>|| I Ibll ||A-7|| + ||V,<I>|| I ||g||J|A'^{I-Pi}/||^ 

< ||V.$||^. \\g\\ ||A-7|| + C ||V.$||^, \\g\\l + ^ ||A-{I - Pi}/||J . 
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For the last term in ()2.94p . we do not need the sphtting and similarly we have 

(A-^ (V,<i> • V,g),A~'f) < \\A-' (|V,<I>||V,<7|2)|| ||A-V|| < l|V.$||^2 HV.^II ||A-7|| • (2.101) 

Consequently, in light of (I2.95p - (l2.10ip and the definitions of Sn and Vjy, we deduce from (|2.94p 
that (l2:92ll holds for s £ (0, 1/2] and that (12:931) holds for s G (1/2,3/2). □ 

2.3. Proof of Theorems ll.il In this subsection, we will combine all the energy estimates that we 
have derived in the previous two subsections and the interpolation between negative and positive 
Sobolev norms to prove Theorem 11.11 Note that for i = 0,1, we can then absorb the right hand 
side of ()2.78p to obtain, by adjusting again the constant in the definition of £i{t). 



d 
lit 



£t+ ||[v'^{i-Pi}/,v^'{i-P2}5 



E<k<N 



e<k<N 



+ 



E 

e<k<N+l 



< for ^ = 0, 1. 



(2.102) 

i+l<k<N 

In particular, taking ^ = in (j2.102p and then multiplying by a large number K > 0, adding the 
resulting estimate with (j2.80p . we obtain 

/ 



d 
di 



K£o + 



\ 



h\+\0\<N 

I/3|>1 



^J{I-P^}f,^2{I-P^}g 



+ CePiv < 0. 



(2.103) 



We may define Cq ^ times the expression under the time differentiation in (j2.103p to be the instant 
energy functional £n, then we have 

^£n+'Dn<0. (2.104) 
dt 

Integrating ()2.104p directly in time, we get ()1.2ip . Hence, if we assume £n{0) < <Jo for a sufficiently 

small ^0 > 0) then a standard continuity argument closes the a priori estimates that <?Ar(t) < 5. 

Thus we can conclude the global solution with the estimate (|1.2ip by the standard continuity 

argument of combining the local existence result and the a prior energy estimates. 

Now we turn to prove (jl.22p and (jl.24p . However, we are not able to prove them for all s at this 

moment. We shall first prove them for s G [0, 1/2]. 



Proof of (ll.22p - (ll.24|) for s G [0, 1/2]. First, integrating in time the estimate (12.921) of LemmaE 
by the bound (|1.2ip . we obtain that for s £ (0, 1/2], 



|A-^/W|| 



< 



|A-Vo||' + C f (||A-^/(r)|| + 1) Vn(.t) dr 
Jo 



1 + sup A" 

0<T<t 



7(r)| 



(2.105) 



This together with (fL2n) gives (fL22D for s G [0, 1/2]. 

Next, we take ^ = 0, 1 in (|2.102p and recall from the definition of the energy functional ££{t) 

that there is only one exceptional term ||V^Pi/(t)|| < ||VV|| that can not be bounded by the 
corresponding dissipation in (I2.102p . The key point is to do the interpolation between the negative 
and positive Sobolev norms by using Lemma lA.4( 

e+a 



< c 



e+i+s 



This together with the bound ()1.22p yields that there exists Cq > such that 



>C||V/|| 



1+ 



e+3 



1+ 



l+a 



(2.106) 



(2.107) 
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Hence, by (|2.1U7I) and (fLUJ, we deduce from (I2.1U2|) that 



—£i + Co{£e] < for £ + s / 0. (2.108) 
at V / 

Solving this inequality directly and (jl.2ip again, we obtain 

£i{t) < (SiiOf^^^^^''^ + Co(£ + s)t) ^^^'^ < Co{l + ty^^+'^ for i + s^O. (2.109) 



Taking £ = 1 in (j2.109p together with ()1.2ip . we obtain (jl.24p . Also note that there is only one 
exceptional term ||Pi/(t)|| in £]\f that can not be bounded by T>n, so using the same arguments 
leading to (|2.108p with i = 0, we can obtain 



£N{t) < (£n{0)-^/' + Costj < Co(l + t)-' for s > 0. (2.110) 
This together with (|1.2ip again gives (|1.23p . □ 
Now we can present the 

Proof of ([L22])-([L2lD for s E (1/2,3/2). Notice that the arguments for the case s £ [0,1/2] can 
not be applied to this case. However, observing that we have /o S L^Hx since L1H~^ n LIL^. C 
LlH-^ for any s' G [0, s], we then deduce from what we have proved for (jl.22p -(|1.24p with s = 1/2 
that the following decay result holds: 

||vV(i)|P < Co(l + t)"(^+^) for ^ = 0,1. (2.111) 

l<k<N 

Hence, by (I2TTB and ([Tn]) . we deduce from (12:931) that for s G (1/2,3/2), 

||A-7(t)||' < ||A-7o||' + C j'^ ((||A-7(r)|| + 1) Vn{t) + ||/(T)f ||V/(r)f-2^) dr. 
<Co + Co sup ||A-V(r)||+Co [\l + tY^'^I^-'Ut 

0<T<t Jo 



<Co{l+ sup A-7(r) 

0<r<t 



(2.112) 



This proves (jl.22p for s G (1/2,3/2), and we may then repeat the arguments leading to (|1.23p - 
(fT:2i]l for s G [0,1/2] to obtain (fOSD - lfOill for s G (1/2,3/2). The proof of Theorem O is 
completed. □ 



3. Weighted energy estimates and proof of Theorem 11.41 

3.1. Weighted energy estimates. In this subsection, we will derive the weighted energy esti- 
mates on the spatial derivatives of the solution, and then with the help of these weighted norms we 
will derive some further energy estimates compared to the basic energy estimates derived in Section 
12.11 The following lemma provides the weighted energy evolution of [{I — Pi}/, {I — P2}5']. 

Lemma 3.1. If £]\f{t) < 5, then for k = 0, . . . , N — 1, we have 

^ II [{I - Pi}/, {I - P2}9] \\l + C\\ HI - Pi}/, ^{I - P2}9] f<V^ (3.1) 
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and 



It 



||[v^{I-Pi}/,V'={I-P2}5 

l<k<N-l 

+ C Yl ||^V'={I-Pi}/,^V^^{I-P2}5 



l<fc<JV-l 



(3.2) 



Proof. We only illustrate the estimate on {I — P2}fl'- Applying d'^ (with < I7I = A; < — 1) to 
(jl.SP n and then taking the inner product with i'd'^{l — P2}5, we obtain 



(3.3) 



+ {dtd^P29 + V • V^d^P29 - 5^V,«> • v^, ud^I - Pals) 
(rig, f) + ^V.c^ • vf - • V./) , i^d^I - P2}9 
By the linear estimate ()A.20|) in Lemma \A.7\ we know 

{€28^1 - P2}9, - P2}9) > \ pd^{l - P2}9f - C W{1 - P2}9\\1 



(3.4) 



(3.5) 



The second line in (13. 3p can be bounded by, as in Lemma 

\ W{1 - P2]9\\l + C \\V.d^9f + c WVM^ ■ 

Now we turn to the third line in (|3.3p . First, by the estimate ()A.23P (with 77 = 0) of Lemma [A. 81 
and applying the similar argument that for Ji, we have 

(a^r(5,/), 1.5^1 -P2}5) = C:i^{^{d^'9,d^-^^f),ud^{l-P2}9) 
< Y ||r(9^^<?,5^-^V)|| \\yd'<{l-P2}g\\ 

71 <7 



71 <7 



< 5 



l<£<Af-l 



(3.6) 



Comparing this with (j2.49p . the only difference is that we replace the z^^/^-weighted norm by the 
I'-weighted norm. This observation is also valid for the estimates of the last two term in the third 
line in (13. 3p . Hence, we may easily complete the estimates for {I — Pa}^ with this repalcement. 
Applying the similar argument and the observation to {I — Pi}/, and summing over I7I = k, we 
get 



2di 



V'={I-Pi}/,V^{I-P2}<7 



1 

+ 4 



< 



+ 



+ 



i/V'={I-Pi}/,z.V^{I-P2}<7 

2 



2<i<N-l 



(3.7) 



+ 6 Y ||[^V^{I-Pi}/,^V^{I-P2}5 



l<^<Af-l 
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Summing the above up k from 1 to — 1, by the definition of Pat, since 5 is small, we obtain (j3.2p . 
The estimate (j3.ip follows similarly and we omit the details. □ 



By ([33]), we know that if £n{0) + ||[/o,5o]||^ and is smah, then S^it) + \\[f{t), g{t)]\\l is small. 
With the help of this weighted bound, we can improve the energy estimates derived in Section [2.11 

Lemma 3.2. If £N{t) + \\[f {t) , 9{t)]\\l < then for k = 0,...,N -1, we have 



d 
It 



+ 



+ c 



V'^{I-Pi}/,V'={I-P2}5 



< 5 



+ 



+ 



2<e<N 



(3.8) 



and for k = N , we have 
d 



dt 



K/, V^] r + II V^V.cI>|n + C II [V^{I - Pi}/, V^{I - F2}9] 



<5(||[V^/,V^]f + Yl ||VV.<I> 

2<e<N 



(3.9) 



Proof. Comparing ()3.8p - ()3.9p with (j2.34p - ()2.35p of Lemma [2.2l the only difference is that we remove 
the last z^^/^-weighted summing term from ()2.34p - (j2.35p . Hence clearly, we only need to improve 
the estimates of those terms in the proof of Lemma 12.21 that leads to this i/^/^-weighted summing 
term. Indeed, these terms are Ji, J2 and /i,/2- We shall only illustrate the improved estimates of 
Ji, and the other terms can be treated in the same way. More precisely, we shall revisit the term 

■■|Vl7ilg|2|V^-|^ll{I-Pl}/|, 



J 



12 



(3.10) 



Note that now we can also bound the i^- weighted factor by the energy, so we can pursue to adjust 
the index. For A; = 0, . . . , — 1, if |7i| =0, then we have 



•^12 < NIl^oz.^ ^^{I-Pi}/ 



< 6 



V^'{I - Pi}/ 



(3.11) 



if I71I ^ li then by Lemma lA.61 and Lemma lA. 21 we have 
^12 < 



< 



















LlLl 





< IIVMI'-^ 



l7l|-l 



71 I - 1 



(3.12) 



[I-Pi}/||. ^ V^I-Pi}/ 



<5 



+ 



V^{I-Pi}/ 



where we have denoted a by 

1 _ |7i| . 
3 3 



l7i| 



'2 3 / k 



fc-(|7i|-l) - 2 



(3.13) 
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Hence, we have that for A; = 0, . . . , — 1, 



J^2 < 5 



+ 



V^I-Pi}/ 



Now for k = N, if |7i| > — 1, by Lemma lA.61 and Lemma \A.2\ we estimate 



Jl2 < 



VlTi 



,l/2yJV-|7i||j_p^|j 



< 



yhilg i.l/2yiV-|7i|{I_p^}J 



2|7l|-3 



< iiv^^r"^^ iiv^5ii ii{i 

<5(||V^5|| + ||V^{I-Pi}/||J 
where we have denoted a by 



I7i| 



axil 



a 



2|7l|-3 1 2|7i|-3 

Pi}/I|. V^{I-Pi}/ l^ 



2|7i|-3A^^^2|,,|-3 



2N 
3N 



2N 



< N: 



2{N-\ji\) + 3 

and if |7i| < — 2, by again Lemma lA.61 and Lemma lA.21 we estimate 



Jl2 < 



vItiI( 



y7V-|7i||j_p^|y < y|7i| 



V^-ItiI{i_Pi}/ 



< iiv-sir^^ iiv^^ii " ii{i-Pi}/ii.~ iiv^^{i-Pi}/ 

<5(||V^5|| + ||V^{I-Pi}/||J, 
where we have denoted a by 



111 ^-I7ll 

AT IIxvjV/t _ -p. \ fll N 



ml 

3 



1 N 
2~ Y 
3N 



N 

3N 



1 a 

2 ~ 3 



X 1 



ml 

N 



a 



IN- < — :— since 7i < — 2. 
2(Ar-|7i|) - 4 - 



Hence, we have that for k = N, 

^i2<5(||V^5|| + ||V^{I-Pi}/||J 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



(3.19) 



Comparing the estimates (j3.14p and (j3.19p with (j2.48p . we have succeeded in removing the 
z/"^/^-weighted summing term from the estimates of J12. Applying the similar argument, we can 
remove this z^^/^-weighted summing term from the estimates of J2 and Ii,l2, and hence we get 



By the estimates p.ip - ()3.2p . we know that if £n{0) + Y2o<k<N-i II [^'^foi '^^do] ||^ is small, then 

£Nit) + Eo<fc<7V-i II [V''/(t), ^^^(t)] \\l is small. With the help of this weighted bound, we can 
deduce a further energy estimate of g and which implies the exponential decay of g and Vx^- 
This energy estimate also can be used to kill the summing term related to Vx^ in the right hand 
side of (|3.8p ~ ()3.9p . and hence the energy estimates in Lemma 13.21 will be improved. 



Lemma 3.3. If SN{t) + '^o<k<N-i II [^'^/C^)) ^'^sl^)] ||^ ^ then there exists an equivalent energy 
functional 

2 



0<k<N~l 



0<k<N 



(3.20) 
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such that the following inequality holds: 

d 



dt 



< 0. 



(3.21) 



^0<A:<Af-l 0<fc<Af 

Proof. The standard V*^ energy estimates on (ll.SP o yields that for k = 0, . . . , N — 1, 



V^{I-P2}5 



< (v'' (r{g, f) + iv.c^ • vf - • V,/) , • 



(3.22) 



To estimate the right hand side of ()3.22p . we will simply bound the norm of the /-related factors 
by 5. More precisely, by the collision estimate (|A.23p (with r/ = 0) and Sobolev's inequality, we 
have 



0<i<k 

^ E \\^'f\l E IK'^ 
E p5 

0<fc<Af-l 

Applying the similar argument to the other two terms, we get 



0<fc<Af-l 
|2 

I v 



V^g 



(3.23) 



d_ 
'dt 



0<k<N-l 



E [\\^'9f + 



^ 0<A:<Af-l 



< 



+ l|V'=P2<?f 



(3.24) 



0<k<N-l 

On the other hand, we may go back to (12.22P to find that for A; = 0, . . . , — 2 
d 



7k+2^ 



dt 



< IIV^I - P2}<7f + l|V^+i{I - P2}5]f + 



(3.25) 



To estimate ||V''9T2 11 IT 

', we simply bound the norm of the /-related terms by 5 to have 



< 



+ liv'-^f 



(3.26) 



0<fc<Af-2 

Then summing up (j3.25p from A; = to A^ — 2, since 6 is small, we have 
I E E l|V'P2<7f+ E ||V'V,.«J>||' < Yl \NHi-P2}9f. (3.27) 

0<k<N-2 0<k<N-l 0<k<N 0<k<N-l 

A suitable linear combination of (j3.24p and (j3.27p gives ()3.2ip . □ 

The following lemma provides the needed estimates for proving the faster decay of the microscopic 
part {I -Pi}/. 
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Lemma 3.4. If £N{t) + 



0<k<N- 



1 II [VV(i), ^=5(0] II' < ^, then fork = 0,...,N-2, 



d 
dt 



V^{I - Pi}/, V^'{I - F^}g] IP + C II [V'^il - Pi}/, V'il - P2}9 

2 



< 



k+l. 



+ 



+ E 



(3.28) 



0<fc<Af-2 

Proof. The proof of p.28p is similar to those of Lemma I3.1H3. 31 and hence is omitted. □ 

3.2. Proof of Theorem 11.41 In this subsection, we will complete the proof of Theorem 1 1.41 by us- 
ing instead the weighted estimates derived in the previous subsection. Recall that all the statements 
of Theorem 11.11 are valid. 

First, letting £n{0) + ||[/o)5o]||^ be sufficiently small, then by the estimate p.ip of Lemma [3Tl 
we have that <SAr(t) + || [/(i), 11,^ is small. Hence, in light of the estimates in Lemma 13.21 we can 
improve the estimates (|2.78p to be that for ^ = 0, . . . , — 1, 

'^-£,+ ^ ||[VHl-Pi}/,VHl-P2}5]||^+ E 

+ ||v''P25 Yl l|v''v^$f<C55 Y IIV^v^^ " 

e<k<N e<k<N+l 2<k<N+l 

Taking ^ = 2 in (j3.29p . we can then absorb the right hand side to obtain, by adjusting again the 
constant in the definition of £2it), 

p2+ Y ||[v'={I-Pi}/,V'={I-P2}5]|[ 

2<k<N ^ 

+ Y llv^Pi/llV Y E ll^''^^^ ^0- 

3<k<N 2<k<N 2<k<N+l 

Then the time differential inequality ()2.108p corresponds to 

Solving this inequality and recalling the definition of £2, we deduce (jl.26p . 

Now, letting £n{0) + Ylo<k<N-i II [^'^fo, ^''do] ||jy be sufficiently small, then by the estimate (j3.2 

of Lemma [3T| we have that £N{t) + J2o<k<N-i II [^'^/(O) ^'^fi'lO] ||^ is small. Then by Lemma 
there exists a constant A > such that 



(3.30) 



(3.31) 



(3.32) 



Solving this inequality and recalling the definition of £g, we obtain (jl.27p . On the other hand, 
adding the estimates ()3.29p and (j3.2ip . since 5 is small, we obtain 

^ + £:,)+ ^ ||vHi-Pi}/ll^+ E llv'Pi/f 

+ E ll^'^^ll + E ii^'^v^^f <o- 



dt 



0<k<N 



0<k<N+l 



Notice that there is only one exceptional term V^Pi/ in ££ + £g that can not be bounded by the 
corresponding dissipation in (j3.33p . But similarly to (j2.108p . we have that for £ = 3, . . . , — 1, 

d 



dt 



£i + £g] + Co[£i + £ 



< 0. 



(3.34) 
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Solving this inequality will in particular give (jl.28p and (|1.29p . 

Finally, applying the Gronwall inequality to (j3.28p . by (|1.27p and (jl.29p . we obtain that for 
k = 0,...,N -2, 



V^{I-Pi}/(t),V^{I-P2}5(i) 



< e 



-ct 



V^'{I - Pi}/o, V^'{I - P2}5o] f + cj^ e-^(*-^) [V^'+V(r), V^'+i5(7 



+ 



0<fc<A'-2 



dT 



< Coe-^* + Co y* e-^(*-^) ((1 + t)-('=+i+'^) + e"^^) dr 

< Coil + 

This in particular gives (|1.3Up . The proof of Theorem 11.41 is completed. 



(3.35) 

□ 



Appendix A. Analytic tools 

A.l. Sobolev type inequalities. We will extensively use the Sobolev interpolation of the Gagliardo- 
Nirenberg inequality. 

Lemma A.l. Let < m, a < i, then we have 

Iiv^/IL, < ||v™/||ii''||vV 

where < 9 < 1 and a satisfy 



L2 



(A.l) 



Proof. This can be found in [15, pp. 125, THEOREM]. 



(A.2) 



□ 



We shall also use the corresponding Sobolev interpolation of the Gagliardo-Nirenberg inequality 
for the functions on x M^. 

Lemma A.2. Let < m, a < i. Let w{v) be any weight function of v, then we have 



where < 6 < 1 and a satisfy 

1 a /I m 



w 



dv 



, (1 - ^) + , 



Proof. For any function f{x,v), by Lemma lA. 11 we have 

l|v"/||^P < ||v-/||i/ vV 



Li 



(A.3) 



(A.4) 



(A.5) 
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Taking the square of (jA.SP and then multiplying by w{v), integrating over M^, by Holder's inequality, 
we obtain 



w 



NV\\hdv< [ w\\V 



26 



LI 



dv 



2(l-6») / 1 

||V"^/||^2 ) (w-^ 
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1-6 



< 



1 

W2 



dv 



(A.6) 



l|v"^/IL. j dvj 

Taking the square root of (|A.6p . we get (|A.3p . 

A. 2. Negative Sobolev norms. We define the operator A^,s G M by 



Ll 



26 



dv 



□ 



(A.7) 



where / is the Fourier transform of /. We define the homogeneous Sobolev space H'^ of all / for 



which ll/ll^fs is finite, where 



|A^/II 



L2 



L2 



(A.8) 



We will use the non-positive index s. For convenience, we will change the index to be "— s" with 
s > 0. We will employ the following special Sobolev interpolation: 



Lemma A. 3. Let s > and l> Q, then we have 



< 



L2 



L2 



A */||^2 ) where 9 



1 



Proof. By the Parseval theorem, the definition of (|A.8p and Holder's inequality, we have 



L2 



< 



L2 



L2 



L2 



1-6 



L2 



L2 



(A.9) 

(A.IO) 

□ 



We shall use the corresponding Sobolev interpolation for the functions on M;^ x MjJ. 
Lemma A.4. Let s > and l> Q, then we have 



|A */|| , where 



'+1 + S 



Proof. It follows by further taking the L^ norm of ()A.9p over Mj^. 



(A.ll) 



□ 



If s G (0,3), A defined by ()A.7p is the Riesz potential. The Hardy-Littlewood-Sobolev 
theorem implies the following L^ inequality for the Riesz potential: 



Lemma A. 5. Let 0<s<3, 1 < p < q < oo, 1/q + s/3 = 1/p, then 

||A-7L, <II/IIl.- 

Proof. See [H pp. 119, Theorem 1]. 



(A.12) 



□ 
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A. 3. Minkowski's inequality. In estimating the nonlinear terms, it is important to use the 
Minkowski's integral inequality to interchange the orders of integration over x and v. 

Lemma A. 6. Let 1 < p < oo. Let f be a measurable function on x M^, then we have 

([ ([ \f{y,z)\dy] dzY < [ ([ \f{y,z)\Pdzy dy. (A.13) 



In particular, for 1 < p < q < oo, we have 



<\\f\\Lm- (A-14) 



T<i TP 

^z^y " " " ^ z 



Proof. The inequality (jA.lSp can be found in [18, pp. 271, A.l], hence it remains to prove (|A.14p . 
For q = oo, we have 



LfLP = sup 



[ \f\^dy] <([ (sup|/|)'d^) =ll/llLgLgc- (A.15) 



For q < oo and hence 1 < q/p < oo, then by ()A.13p . we have 



tQ rP 

l^z^y 



□ 

A. 4. Boltzmann collision operators. First, we collect some useful estimates of the linear colli- 
sion operators. 

Lemma A. 7. For i = 1,2, we have 

{dhi, h2) = {hi, CM), {Cih, h) > 0, (A.17) 
Cih = if and only if h = F^h. (A. 18) 

Moreover, there exist (Tq > and constants C, C|^| > such that 

{Cih,h)>ao\{I-Pi}h\l, (A.19) 
{udpdh, dph) > I \udph\l - C\h\l, (A.20) 



2 

{dpCih,dph)>\\dph\l-C\p\\h\l. (A.21 



Proof. We refer to ^ Lemma 1] or [22t Lemma 3.1] for (jA.17|) - (|A.18j) . While for the proof of 
([A:T9]) - ([A:2TT) . we refer to [H Lemma 3.2-3.3]. □ 

Now, we collect some useful estimates of the nonlinear collision operator. 

Lemma A. 8. There exists C > such that 

Kr(/li,/l2),/l3)| + |(r(/l2,/il),/i3)l <Csup{l/3/i3}|/j,|2|/i2|2. (A.22) 

V 

Moreover, for any < < 1, we have 

|l/-''r(/ll,/l2)|2 <C{\u^~'^hi\2\h2\2 + \v^-'^h2\2\hl\2] ■ (A.23) 

Proof We refer to O Lemma 2.3] for (1X22]) . and [211 Lemma 2.7] for (IXMj) . □ 
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